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microlensing parallax, orbital motion and detection of blended light 
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ABSTRACT 



Microlensing detections of cool planets are important for the construction of an 
unbiased sample to estimate the frequency of planets beyond the snow line, which is 
where giant planets are thought to form according to the core accretion theory of planet 
formation. In this paper, we report the discovery of a giant planet detected from the 
analysis of the light curve of a high-magnification microlensing event MOA 2010-BLG- 
477. The measured planet-star mass ratio is q = (2.181 ±0.004) x 10~ 3 and the projected 
separation is s = 1.1228 ± 0.0006 in units of the Einstein radius. The angular Einstein 
radius is unusually large #e = 1-38 ± 0.11 mas. Combining this measurement with 
constraints on the "microlens parallax" and the lens flux, we can only limit the host 
mass to the range 0.13 < M/M Q < 1.0. In this particular case, the strong degeneracy 
between microlensing parallax and planet orbital motion prevents us from measuring 
more accurate host and planet masses. However, we find that adding Bayesian priors 
from two effects (Galactic model and Keplerian orbit) each independently favors the 
upper end of this mass range, yielding star and planet masses of M* = 0.671q'^3 M & 
and m p = 1.5lg'| -^JUP at a distance of D = 2.3 ± 0.6 kpc, and with a semi-major 
axis of a = 2±l AU. Finally, we show that the lens mass can be determined from future 
high-resolution near-IR adaptive optics observations independently from two effects, 
photometric and astrometric. 

Subject headings: gravitational lensing: micro - planetary systems 
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Introduction 



Gravitational microlensing i s an important method to detect extrasolar planets (IMao Paczvhski 



1991 



Bennettll2008l ; lGaudill201Cll ). The method is sensitive to planets not eas ily accessible to other 



methods, in particu lar cool and small planets at or beyond the snow line (jBeaulieu et al 



Bennett et al.1 12008 ) , and free-floating planets (jSumi et al.1 12011 



locat i on in the protoplanetary d isk beyond which ices can exist ( 



2006 



The snow line represents the 



2007 



Lecar et al.1 120061 : iKe nnedv et al 



Kennedy k, Kenvonl 12008) and thus the surface density of solids is hi ghest (|Lissauerill987l ). 



According to the core accretion theory of planet formation (ILissaueii Il993l ). the snow line plays 
a crucial role because giant planets are thought to form in the region immediately beyond the 
snow line. Therefore, microlensing planets can provide important constraints o n planet formation 
theories, in part i cular by measuring the mass function beyond the snow line (jGould et al.l 120101 ; 
Sumi et alJlioiol : ICassan et allboij i. 



A major component of current planetary microlensing experiments is being carried out in 
survey and follow-up mode, where survey experiments are conducted in order to maximize the event 
rate by monitoring a large area of the sky one or several times per night, while follow-up experiments 
are focused on events alerted by survey observations to densely cover planet-induced perturbations. 
In this mode, high-magnification events are important targets for follow-up observations. This is 
because the source trajectories of these events always p ass close to the central p e rturbation region 
and thus the sensitivity to planets is extremely h igh (jGriest Safizadehl Il998l ; iRhie et al.l 120001 ; 
Rattenburv et alj 120021 : lAbe et al.ll2004l ; lHanll2009l ). In addition, the time of the perturbation can 
be predicted in advance so that intensive follow-up observation can be prepared. This leads to an 
observational strategy of monitoring high-magnification events as intensively as possible, regardless 
of whether or not they show evidence of planets. As a result, the strategy allows one to con struct 
an unbiased sample to derive the frequency of planets beyond snow li ne (Gould et al.ll2010f l. For 
the alternative low-magnification channel of detection, see for instance ISumi et al.1 (|2010i ). 



In this paper, we report the discovery of a giant planet detected from the analysis of the light 
curve of a high-magnification microlensing event MOA 2010-BLG-477. Due to the high magnifica- 
tion of the event, the perturbation was very densely covered, enabling us to place constraints on the 
physical parameters from the higher-order effects in the lensing light curve induced by finite source 
effects as well as the orbital motion of both the lens and the Earth. We provide the most probable 
physical parameters of the planetary system, corresponding to a Jupiter-mass planet orbiting a 
K dwarf at about 2 AU, the system lying at about 2 kpc from Earth. 



Observations 



The event MOA 2010-BLG-477 at coordinates (RA, DEC)=(18 h 06 m 07.44 s , -31°27'16.12") 
(J2000.0), (l,b) = (0.046°, —5.095°) was detected and announced as a m i crolensing alert e vent 



by the Microlensing Observations in Astrophysics (|MOA: Bond et al. 



2001 



Sumi et al J 12003 ) col- 
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laboration on 2010 August 2 (HJD' = HJD - 2450000 = 5410.9) using the 1.8 m telescope of 
Mt. John Observator y in New Zealand. Th e event was also observed by the Optical Gravitational 
Lensing Experiment (jOGLE: Udalskil 120031 ) using the 1.3m Warsaw telescope of Las Campanas 
Observatory in Chile. 

Real-time modeling based on the rising part of the light curve indicated that the event would 
reach very high magnification and enabled us to predict the time of the peak. Followed by this sec- 
ond alert, the peak of the light cur ve was densely covered by follow-up groups of the Probing Lensing 
Anomalies NETwo rk (PLANE T) (beaulieu et alJbood ). Microlensing Follow-Up Network (//FUN) 
(jGould et al.l 1200a ) . RoboNet (jTsapras et al.l 120091) . and Microlens ing Network for the Detection 
of Small Terrestrial Exoplanets (MiNDSTEp) (|Dominik et al.ll2010l ). More than twenty telescopes 
were used for the follow-up observations, including PLANET 1.0m of South African Astronomical 
Observatory (SAAO) in South Africa, PLANET 1.0m of Mt. Canopus Observatory in Tasmania, 
Australia, PLANET 0.6m Perth Observatory in Australia, PLANET 0.4m ASTEP telescope at 
Dome C, Antarctica, /iFUN 1.3m SMARTS telescope of the Cerro Tololo Inter-American Observa- 
tory (CTIO), Chile, /iFUN 0.4m of Auckland Observatory, /iFUN 0.36m of Farm Cove Observatory 
(FCO), /iFUN 0.36m of Kumeu Observatory, /iFUN 0.4m of Possum Observatory, /iFUN 0.4m of 
Vintage Lane Observatory (VLO), /iFUN 0.3m of Molehill Astronomical Observatory (MAO), all 
in New Zealand, /iFUN 0.46m of Wise Observatory, Israel, /iFUN 0.8m of Teide Observatory at 
Canary Islands (IAC), Spain, /iFUN 0.6m of Pico dos Dias Observatory, Brazil, RoboNet 2.0m 
of Faulkes North (FTN) in Hawaii, U.S.A., RoboNet 2.0m of Faulkes South (FTS), Australia, 
RoboNet 2.0m of Liverpool (LT) at Canary Islands, Spain, MiNDSTEp 1.54 m Danish telescope at 
La Silla, Chile, and MiNDSTEp 1.2m Monet North telescope at McDonald Observatory, U.S.A.. 

Please note that this paper reports the first microlensing observations from Antarctica. Un- 
fortunately, the data quality of these pioneering observations was not high enough to be included 
into the models, but we give a short overview in Appendix [AJ In order to better constrain the 
second-order effects, new observations were taken at the /iFUN 1.3 m SMARTS telescope at CTIO 
during the 2011 campaign. 

To better characterize the lensed source star, spectroscopic observations were conducted near 
the peak of the event (HJD' = 5422.5) by using the B&C spectrograph on the 2.5 m du Pont 
telescope at Las Campanas Observatory in Chile. The resolution was R = 1400, corresponding to 
AA = 4.6 A. A comparison between the observed spectrum and synthetic ones was conducted to 
derive the effective temperature, gravity and metallicity of the source star. The usable part of the 
spectrum is only ~ 1000 A due to some scattered light issues with the instrument, particularly a 
problem for this faint star. We started by fitting standard stars using Ha, MgB, and NaD, but 
found that fits with Ha alone resulted in the most accurate temperatures, so this is the diagnostic 
we used. The derived parameters of the source star are T e g- = 5950 ± 150 K for log g = 4.0 and a 
solar metallicity. The corresponding fit to the Ha line is shown in Figure HJ 



The data collected by the individual groups were initially reduced using various photometry 
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codes developed by the individual groups. For the data sets from SAAO, FTS, Possum, Canopus, 
Perth, Danish and Monet North, we use photometry from rereductions obtained with the pySIS 
package, described in more details in Appendix [Bl 

Figure [2] shows the light curve of the event. The event is highly magnified with a peak 
magnification A ~ 400. Outside of the region H JD' = [5417,5425], the light curve is consistent 
with a standard single-lens curve (jPacz vriskil 1 1 9861 ) . The perturbation is composed of two spikes at 
HJD' ~ 5420.4 and 5420.9 and two bumps at HJD' ~ 5421.0 and 5422.4. 



3. Strategic Overview 

As in the great majority of planetary microlensing events, we are able to measure the "angular 
Einstein radius" (projected on the sky) #e ; but not the projected Einstein radius (projected on the 
observer plane) f E - Consequently, it is challeng ing to estimat e the physical parameters of the lens. 



The two radius quantities are related to those (|Gouldll2000bl . see) by 



e E = K7T m iM; vr E = — — ; r E = (1) 

K M 7TE 

where k = 4GM/AUc 2 ~ 8.1439 mas Mq 1 , M is the lens mass in M , and 7r re i = AXJ(l/D L -l/Ds) 
is the lens-source relative parallax (Dl an d Dg are the lens and source distances). Since #e is well 
measured, the product 7r re i M is also well determined, but in the present case, the ratio 7r re i/M 
is poorly constrained, and so it is difficult to estimate M alone. In this section, we provide an 
overview of the various techniques that we use to place constraints on the individual quantities 7r re i 
and M. 

We will show that here #e is large enough to enable a substantial constraint from upper limits 
on the lens flux. That is, from Equation ((!]), large # E implies large M or 7r re i, both of which lead to 
brighter lenses. This will lead to the unambiguous conclusion that the lens is nearby, Dl < 3kpc, 
with mass M < 1 Mq. 

The light curve of this event enables stronger constraints than is usually the case because we are 
able to obtain a measurement of one of the microlens parallax component. The microlens parallax 
is actually a vector, 7T E ; with the magnit ude given by Equation ([TJ and the direction by the lens- 



source relative motion (jGould et al.lll994l ). One component of the parallax only weakly constrains 
the scalar 7Te, but it constrains the direction of proper motion, which will be very important for 
future observations (see below). Moreover, one parallax component actually does provide a robust 
lower limit on the mass M > 0.13 Mq and distance Dl > 0.5 kpc. 

To proceed further, we must apply a Bayesian analysis, but here as elsewhere we are more 
fortunate than is typical. As always, there are Bayesian priors from a Galactic model, and in this 
case these strongly prefer the upper end of the mass/distance range permitted by the light curve. 



But in addition, the reason that only one dimension of 7Te can be robustly measured is that the 
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other dimension is degenerate with orbital motion of the planet (jBatista et al.ll201ll ). Thus, it is 
necessary to fit simultaneously for both parallax and orbital motion. Bayesian priors on the orbital 
motion then independently also prefer the upper end of the allowed mass/distance parameter space. 

Finally, we predict that high-resolution imaging could measure the mass and distance of the 
lens through two independent effects, photometry and astrometry. The first is obvious: different 
mass/distance combinations yield different color /magnitude measurement. The main point here 
is that the large value of #e virtually guarantees that the lens will be detectable. The second is 
less obvious: different mass/distance combinations also predict different directions of source- lens 
relative proper motion, which can be measured as the lens and source are separated over the next 
several years. 



4. Modeling 

4.1. Treatment of photometric errors 

The photometric error bars of the various data sets from individual observatories are generally 
not accurate enough to be taken at face value. They need a rescaling to reproduce the dispersion 
of contiguous data points in a given night. At the same time, outliers must be identified and 
removed. These are important preliminary steps to modeling, because the weight of a given data 
set to constrain the model depends on how large its error bars are, compared to other data sets. We 
describe our adopted noise model and rescaling factors for each observatory data set in Appendix [Cl 



4.2. Static Binary Model 

We first test a static binary lens model. The corresponding parameter set includes the three 
single lens parameters: the Einstein time scale, £e, the time of the closest lens-source approach, to, 
and the lens-source separation at that moment, uq, and it also includes the three binary parameters: 
the mass ratio of the companion to its host star, q, the projected separation between the lens 
components in units of the Einstein radius, s, and the angle between the source trajectory and 



the modification of magnifications caused by the finite-source effect ( 


Nemiroff & Wickramasinghe 


1994: 


Witt & Mao 


1994: 


Gould 


1994: 


Bennett & Rhie 


1996: 


Vermaak 


2000). This requires us to 



include an additional parameter: the normalized source radius 



P* = 7T 



0* 
0e 



(2) 



where 6* is the angular source radius. Evaluation of p* from the model together with the measure- 
ment of 6* (see Section ETTj) will then yield 6e and thus the product n re \M (see Equation [1]). 
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4.3. Finite-source Effect 



Since the event MOA 2010-BLG-477 involves caustic crossings and approaches, one must com- 
pute magnification affected by the finite-source effect. Computation of finite magnifications is based 
on the numerical ray-shooting method. In this method, a large number of rays are uniformly shot 
from the image plane, bent according to the le ns equation, and land on the source plane. The lens 
equation is represented by (|Witt fc Madll995l ) 

" m k /M 



c 



k=l 



ZL,k 



(3) 



where £ = £ + irj, Zl^ = Xlu + WLki an d z = x + iy are the complex notations of the source, 
lens, and image positions, respectively, z denotes the complex conjugate of z, and are the 
masses of the individual lens components, and M is the total mass of the lens system. Then, the 
magnification is computed as the ratio of the number density of rays on the source plane to the 
density of the image plane. For the initial search for solutions in t he space of the gr id parameters, we 
accelerate the computation by using the "map making" method (IDong et al.ll2006l ). In this method, 
a magnification map is made for a given set of (s, q) and then it is used to produce numerous light 
curves resulting from different source trajectories instead of re-shooting rays all over again. We 
further accelerate the comp utation by using a semi-analytic hexade capole approximation for finite- 



magnification computation (|Pejcha &: Heyrovsky 



2009 



Gouldll2008l ) in the region where the source 



location is not very close to the caustic. In computing finite magnifications, we consi der the effect 
of limb darken ing of the source star surface by modeling the specific intensity as (jMilne 
An et alJliooa ) 



1921 



T - Fx 



1 COS ( 

2 



(4) 



where T\ is a limb-darkening coefficient (hereafter LDC), Fx is the total flux from the source star, 
cf) is the angle between the direction toward the observer and the normal to the stellar surface. 
From the x 2 improvement we find that the limb-darkening effect is clearly detected. We compute 
the LDCs for Equation 0] as accurately as possible, including a proper treatment of the effect of 
extinction. We refer the interested reader to Appendix [D] for details. 



4.4. Microlensing Parallax and Planet Orbital Motion 



We then test if second-order effects are present in the residuals of t he light cur y e. These effects 
may have several origins: orbital motion of the Earth around the Sun (|Gouldlll992l ). which induces 
a deviation of the lens-source motion from rectilinear, orbital motion of the planet about the lens 
star, orbital motion of the source star if it is a binary, and possible additional objects (planets or 
stars) in the lens system. 

In fact we will limit our study to the first two effects, namely microlens parallax and planet 
orbital motion. A binary source in which the companion is either not lensed or too faint to 
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contribute to the light curve can mimic the parallax effect due to the acceleration of the sour ce, in 
a so-calle d "xallarap" e f fect. This effect is described for instance in iPoindexter et al.l (|2005l ) and 
tested in iMivake et al.l (|2012l ). However, it implies in its simplest form 5 additional parameters 
versus 2 each for parallax and planet orbital motion. Given the fact that our detection of second- 
order effects, although clear, remains relatively marginal, it would be difficult to trust 5 additional 
parameters constrained by the light curve residuals from a standard model. Moreover, binary stars 
generally have long periods, the effect of which will remain undetected in the short lapse of time 
when our events are observed (months versus years). 

Testing for a third body, either from a circumbinary planet or a second planet orbiting the 
lens star, is beyond the scope of this paper, and again would imply a large number of additional 
parameters to fit, without any guarantee to obtain mea ningful results. In the short history o f 



planetary microlensing, it has been evidenced only once (jGaudi et al 



2008 



Bennett et all 1201(3 ) 



and in that case only because a feature in the light curve could not be fitted by a single planet model. 
There is no such feature in the present light curve. As the domain of possible configurations to test 
if we add a third body is vast, the effort expended would be incommensurate with the potential 
scientific return, at least until we get independent evidence that such a modeling effort is necessary: 
see Section [6j 

As discussed in Section [3l the microlens parallax is characterized by a 2-dimensional vector 7Te 
whose magnitude is given by Equation (pQ) and whose direction is that of the lens-source relative 
motion on the plane of the sky. Hence, there are two parameters 7Te = (7Te,n ; tte,e)) the components 
of this vector in equatorial coordinates. 

The planet orbital motion affects the light curve in two different ways. First, it c auses the 
binary axis to rotate or, equivalently, makes the source tra jectory angle change in time (jDominik 



1998 



Ioka et al. 



1999 



Albrow et al. 



2000 



Penny et al.ll201ll ). Second, it causes the projected binary 



separation, and thus the magnification pattern, to change in time. Considering that the time scale 
of the lensing event is of order a month while the orbital period of typical microlensing planets is 
of order years, the rates of change in a and s can be approximated as uniform during the event 
and thus the orbital effect is parameterized by 



a(t) = a(t rei ) +lo x (t — i re f ) 



(5) 



and 



s(t) = s(i re f) + s x (t - t ie {) 



(6) 



where u and s are the rates of change in the source trajectory angle and pr ojected binary separ ation 



in units of yr 1 , respectively, and i re f is a reference time. As explained in lBatista et al. ( 2011 ). the 



effect of planet orbital motion is similar to the microlensing parallax effect in the sense that the 
deviations caused by both effects are smooth and long lasting. Then, if the deviation caused by the 
planet orbital motion is modeled by the parallax effect alone, the measured parallax would differ 
from the correct value. Therefore, the orbital motion effect is important not only to constrain the 
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orbital properties of the lens, but also to precisely constrain the lens parallax and thus the physical 
parameters of the lens system. 

For each tested model, we search for the solution of the best-fit parameters by minimizing 
X 2 in the parameter space. In the initial search for solutions, we divide the parameters into two 
categories. Parameters in the first category are held fixed during the fitting, and parameter space 
is searched with a grid. For the parameters in the second category, solutions are searched by using 
a downhill approach. We choose s, q, and a as the grid parameters because they are related to the 
light curve features in a complex way, where a small change in the parameter can result in dramatic 
changes in the light curve. On the other hand, the remaining parameters are more directly related 
to the features of the light curve and can thus be searched by using a downhill approach. For 
the downhill \ 2 minimization, we use a Markov Chain Monte Carlo method. Brute-force search 
throughout the grid-parameter space is also needed to test the possibility of the existence of local 
minima that result in degenerate solutions. For the light curve of MOA 2010-BLG-477, we find 
that the other local \ 2 minima have \ 2 values much larger than the best fit by Ax 2 ~ 6000 and 
are therefore not viable solutions. Once the space of the grid parameters around the solution is 
sufficiently narrowed down, we allow the grid parameters to vary in order to pin down the exact 
location of the solution and to estimate the uncertainties of the parameters. 



Modeling was also done independently using the method of iBennettl (120101 ). and this analysis 



reached the same conclusions. This independent analysis also uses a slightly different implementa- 



tion of the planetary orbital motion parameters (IBennett et al.ll2010l ). Equations [5] and [6] describe 



the orbital velocities, which are the first order contribution of orbital motion. To second order, 
we have only one component of acceleration because this must be directed towards the host star. 
But, one additional parameter is also all that is needed to describe a circular orbit, so we can 
add the planetary orbital period, parameter and replace the constant velocities (in polar 

coordinates) of Equations [5] and [6] with the projection into the plane of the sky of th e circular orbit 



descr ibed by s, a, s, u, and T or b. Unlike the case of OGLE 2006-BLG-109Lb,c (jBennett et al 



20101 ). the value of T or b does not have an influence on the light curve model x 2 values. 



However, T or b is still useful because it can be used to help constrain the other orbital parameters 
to values consistent with a physical orbit. (This is an issue because it is quite possible to have s 
and to values that are not consistent with a bound orbit, and this is a simple way to ensure that 
this is not the case). If we assume that is known, then we can calculate the lens system mass, 
M = #e/(k7Te), which follows from Equation [TJ This also allows us to determine 7r re i, but we 
cannot determine the lens and source distances, Dl and D$, separately. However, we can use 
the T or b value to determine the orbital semi-major axis, under the assumption that the orbit is 
circular. This allows us to determine the Einstein radius, Re = DlOe in physical units, and since 
R\ = AGMD\D git veil '{c 2 : AU), we have a second relation between D$ and Dl- But, since the 
source star is very likely to be in the Galactic Bulge, we also have approximate knowledge of D$- 
Therefore, we can apply a constraint on the value of D$ implied by the light curve parameters, 
D s = 8.0 ± 1.2 kpc. 
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Table 1. Best-fit Model Parameters 



parameter 


standard 


parallax 


model 
orbit+parallax (uq < 0) 


orbit+parallax (un > 0) 


x 2 


6420.051 


6411.450 


6369.345 


6365.664 


t (HJD') 


5420.93685 


5420.93702 


5420.93773 


5420.93915 


u 


-0.003562 


-0.003540 


-0.003372 


+0.003404 


t E (days) 


42.55 


42.94 


46.57 


46.92 


s 


1.12282 


1.12338 


1.12385 


1.12279 


q 


0.0023808 


0.0023694 


0.0022075 


0.0021809 


a (rad) 


3.67605 


3.67528 


3.68023 


2.60087 




0.0006429 


0.0006403 


0.0005861 


0.0005764 


fs 


0.4111 


0.4082 


0.3757 


0.3731 




997.3 


997.9 


1045.8 


1059.8 






+0.37 


+0.27 


+0.77 






+0.012 


+0.02 


-0.11 


s (yr^ 1 ) 






+0.64 


+0.86 


w (yr" 1 ) 






+0.06 


-1.28 


^orb (yr) 






14.25 


4.59 


tref (HJD') 




5421 


5421 


5421 



Note. — HJD' = HJD — 2450000. t re f is the reference time of the model, when the model 
reference frame moves at the same speed as the Earth and a(t) = a and s(t) = s. 
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For this event, T or b is not really constrained by the light curve measurements, so this constraint 
serves to force T or b toward a value consistent with a circular orbit for a Bulge source. This constraint 
on the source distance also serves to enforce a constraint on the orbital velocity parameters, s and 
u>. They must also be consistent with a circular orbit for a Bulge source. Parameters that satisfy 
this constraint are also consistent with most orbits with moderate eccentricity, e ;$ 0.5. But, these 
parameters are not consistent with orbits with the highest possible transverse velocities. In fact, 
the light curve measurements marginally favor implausibly large s and u values corresponding to 
orbits that are either unbound or just barely bound. These barely bound orbits with large s and 
uj values have high eccentricities that just happen to have been observed with motion in the plane 
of the sky during the brief time near periapsis. The best fit model that is consistent with a bound 
orbit is such a model, which has an orbit with an extremely low a priori probability. This low a 
priori probability makes such a model much more unlikely than the best fit model with the circular 
orbit constraint, so we report the best fit model with the circular orbit constraint as the "best fit 
model" in Table [TJ 

However, the unlikely models with values of s and u) do tend to have better light curve x 2 
values, and they also cover a large volume of parameter space. So, they should not be ignored in our 
consideration in the range of possible physical parameters for the MOA 2010-BLG-477L planetary 
system. Therefore, we do not enforce this source distance constraint in the MCMC runs that we 
use to estimate the distribution of likely physical parameters for this system. Instead, we apply a 
more general constraint on the orbital and Galactic parameters of the lens system as discussed in 
Appendix El 



5. Results 

In TableHJ we present the results of modeling along with the best-fit parameters for the 3 tested 
models. The best-fit light curve is presented in Figure [2j In Figure El we also present the geometry 
of the lens system. It is found that the perturbation near the peak of the light curve was caused by 
the source crossings and approaches of the caustic produced by the binary system with a low-mass 
companion. The measured mass ratio between the lens components is q = (2.181 ± 0.004) x 10 -3 
and thus the companion is very likely to be a planet. The measured projected separation between 
the lens components is s = 1.1228 ± 0.0006, which is close to the Einstein radius. As a result, 
the caustic is resonant, implying that the caustic forms a single closed curve composed of 6 cusps. 
The perturbations at HJD' ~ 5420.4 and 5420.9 are produced by the source crossings of one of the 
star-side tips of the caustic. The bumps at HJD' ~ 5421.0 and 5422.4 are caused by the source's 
approach close to the weak and strong cusps on the side of the host star, respectively. 
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5.1. Constraints from Measurement of #e 

The most important constraints on the lens mass and distance come from the measurement of 
#Ej which can be rewritten from Equation ([2]) as 

fe _ «. _ VfJWD 

P* P* 

where Fj is the intrinsic (dereddened) flux of the source and Sj is its intrinsic /-band surface 
brightness. A key point is that the surface brightness does not depend at all on the microlens 
model (just on the source color and/or spectrum). Hence, #e depends on the model only through 
the parameter combination yffl/ P*-, where f s is the instrumental source flux (/ s = 1 corresponds 
to magnitude 18). It is clear from Table Q] that this parameter combination varies very little. 
The biggest uncertainties are therefore in measuring the surface brightness and measuring the 
offset between the instrumental flux f s and the dereddened source flux Fj. Traditionally, these 
are measured simultaneously by determining the offsets of the color and magnitude, respectively, 



of the source from the clump on an instrumental CMD (lYoo et al.ll2004l ). An instrumental CMD 
from CTIO is shown in Figure HI From it, we can read instrumental magnitude in I of 19.07. The 
instrumental V — I color is measured more accurately from a regression of V flux vs. I flux, and 
gives —0.35. 

Comparing to the instrumental clump position, we find (V — I)o = 0.55 ±0.05 and Io = 17.61 ± 
0.15, assuming [(V — I), I]o,clump = (1-06, 14 .42). The color erro r is determined empirically from a 



sample of microlensed dwarfs with spectra (jBensby et al.ll201ll ). while the magnitude error comes 



primarily from the error in fitting the clump centroid and the assumed Galactocentric distance of 
8 kpc (both about 0.1 mag). These lead to an estimate 9* = 0.79 db 0.06. 

We can compare the dereddened color estimate to the color deduced from the high-resolution 
spectrum reported in Section [2j from which we measured an effective temperature of T = 5950 ± 
150K, which corresponds to {V— J)o = 0.65±0.04. These determinations are marginally consistent. 
We adopt the first to maintain the general practice of microlensing papers, but note that if we 
adopted the mean of the two estimates, the inferred value of 0* (and so #e) would rise by 2%. 

Finally, we evaluate 

E = 1.38 ±0.11 mas (8) 

This value is unusually large and implies that the lens must be very massive or very close. Specifi- 
cally 

M vr re i _ 0| 



Mq mas k Mq mas 



0.233 ± 0.036 (9) 



5.2. Constraints from Lens Flux Limits 



The model gives a measurement of the light coming not only from the source, but also from 
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any other stars in the same point spread function, generally called blended light. This light may 
come from the lens itself, a companion to the lens, a companion to the source, or any unrelated 
star on the same line of sight, but not participating to the amplification process. 

The OGLE-III image of this field, displayed in Figure El shows two stars at the target position. 
Their /-band magnitudes in the OGLE-III photometric catalogue is / = 17.446 ± 0.052 for the 
brighter star, numbered 119416, and I = 17.807 ± 0.059 for the fainter one with number 119534. 
They are separated by 1.26", enough to be separated by PSF photometry at good sites such as 
CTIO or Las Campanas (OGLE) in Chile. A difference image analysis (DIA) made on CTIO 
images shows that the microlensed star does not correspond to the position of either star (red circle 
in Figure [5]), but is displaced by 1 CTIO pixel (0.37") from the brighter star. 

Now, the blended light in /-band as measured from DoPhot CTIO photometry by the model is 
Jb = 17.443 ±0.031, and this corresponds precisely to the flux from the brighter star among the two 
OGLE stars, which is not separated from the microlensed target at the scale of the CTIO seeing 
(typically 1"). Knowing that the blended flux comes from an unrelated blended star, the light from 
the lens must be smaller. We can rigorously conclude that the lens has less than half the light in 
the observed blend. Otherwise, the lens and blend would be separated by at least 2 x 0.37" = 0.74" 
in the OGLE image, and so would have been at least marginally resolved. 

An additional argument showing that the lens is faint enough to remain undetected comes from 
its large relative proper motion (10.3 ±0.8 mas/yr). DIA analysis of OGLE-III good seeing images 
separated by 3.3 years shows no residual at the target position, which proves that no detected star 
has moved during this period. 

Combining this limit with Equation (J9j) yields strong constraints on the lens. The lens must 
be closer than the source, and so be at or closer than the Galactocentric distance, and suffer the 
same or less extinction. The se imply It ;i o = Ih — 2.51og(0.5) — Ai > 17.0, and so Mj > 2.5, which 



corresponds to M < 1.4 M© (jStraizvs Kuriliendll98ll ; iBessell Brettlll988l ) . Then, even the 2 a 



limit from Equation ([9]) implies 7r re i > 0.10 mas and so Dl < 4.4kpc. But Jl,o > 17.0 then implies 
Mj > 3.8, which corresponds to M < 1.0 M©. Cycling through this argument one more time yields 
Dl < 3.7kpc, Mj < 4.2, M < 1 M©. If the lens were in front of some of the dust, this argument 
would become still stronger. However, at the relatively high latitude of this field (b = —5°), most 
of the dust probably lies in front of 3 kpc, and in any case, there is no basis for adopting a more 
optimistic assumption about the dust. 

5.3. Constraints from the Microlens Parallax tie 

Table Q] shows that parallax alone improves the fit by A% 2 = 8.6, while including both parallax 
and lens orbital motion improves it by A% 2 = 54.4. However, while detection of these effects is 
therefore unambiguous, we cannot fully disentangle one from the other. For each effect, one of its 
two parameters is well determined w hile the other i s high ly degenerate with one parameter from 



the other effect. As first discovered by lBatista et all (|201ll ) and further analyzed by lSkowron et al 
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(|201 ll ) . weak detections of parallax and orbital motion lead to a strong degeneracy between 7Tej_ 
(the component of 7Te perpendicular to the projected position of the Sun) and oj (the component 
of orbital motion perpendicular to planet-star axis), while tteji and s are well constrained. The 
impact of this is illustrated in Panel (a) of Figure [6j which shows \ 2 values by color as a function 
of (tte.NjTTe.e)- 

First, it is clear that the parallax vector is almost completely degenerate along a line that is 
8.4° west of north, within 1° of the predicted orientation of tte,± (7.7°). Second, the light curve 
excludes 7Te > 1.3 at 3a. As we have 

M 6v 

— x7r E = — — = 0.169 ±0.013 (10) 

M KM mas v ; 

this corresponds to M > 0.13 M®. Thus, combining constraints from this section and from Sec- 
tion 15. 2\ and using Equation ([9]) , we have 

O.13M <M < l.OM ;O.5kpc < D L < 2.8kpc (11) 

Note that, because the parallax contours pass through the origin, parallax provides no additional 
constraint at low 7Te, i.e., at high mass. 

A similar diagram is given for the slightly disfavored uq < solution in Figure [71 



5.4. Post-Bayesian Analysis 

Equation (|lip defines the limits of what can be said rigorously about the lens mass and distance 
based on current data. As we discuss in Section [6l the lens could almost certainly be detected by 
high-resolution imaging, which would probably completely resolve the uncertainty in Equation (|lip . 
In the meantime, we can perform a Bayesian analysis based on a Galactic model and constraints 
from a Keplerian orbit. Figure [6] shows that each of these constraints separately favors small values 
of 7Te and hence, relatively large masses and distances (within the limits set by Equation |11|). 
and so tend to reinforce each other. Panel (a) shows the raw results of the MCMC. Panel (b) 
shows the same chain, post-weighted by the Bayesian prior due to the Galactic model and the 
flux constraint. The latter, which implies M < 1.0 M (Section \5.2\i is r esponsible for t he cir cular 



"holes" at the centers of panels (b) and (d) . As discussed in detail by iBatista et al.1 (|201ll ) this 
includes terms not only reflecting the density of lenses along the sight and the expected distribution 
of proper motions, but also a Jacobian transforming from the "natural microlensing variables" to 
the physical system of the Galaxy. It is this last term that actually dominates, particularly in 



the "n ear field" {Dl < 2.8 kpc) permitted by the constraints, and is roughly ~ D\ (jBatista et al 



(|201ll ). Equation [18]). This weighting is so severe that one must ask whether the result of the 
weighting is plausibly compatible with the raw \ 2 from the light curve. In this case, the solutions 
favored by the Bayesian post-weighting are compatible with the raw x 2 minimum at better than 
2<7, so there is no major conflict. 
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Panel (c) shows t he results of post-we ighting the chain only by the orbital Jacobian. This is 



described in detail bv lSkowron et al.l (|201ll ) for the case of chains with complete orbital solutions. 
These contain two orbital parameters (called 7 2 and s z ) in addition to the two first-order orbital 
parameters considered here (s and oS). These higher order parameters would be completely uncon- 
strained in the present problem, so we simply resample the chain with a uniform integration over 
these two parameters. Panel (c) clearly also favors more massive, more distant lenses (although not 
the most massive), but it is not immediately obvious why. Appendix lEl details the reasons behind 
this result. Finally, we note that, as expected, the combined effect of these two priors shown in 
Panel (d) is stronger than either separately. 

Why do the Galactic model and Kepler constraints each favor more distant (and more massive 
lenses) than the light curve alone (see Fig. [6])? The Galactic model constraint is virtually guaranteed 
to favor more distant lenses because, as mentioned above, most of the weighting is simply due to a 
coordinate transformation from microlensing to physical coordinates, and there is more phase space 
at larger distances. Since the distance errors are fairly large, this effect will be relatively strong. 
By contrast, the Kepler constraint could have just as easily favored more-distant as less-distant 
lenses. The most likely explanation is then, simply, that the light curve prediction of the distance 
is too close by 1.7 sigma, and the constrained value (panel d) is a better estimate. This conjecture 
is testable by future observations, as described in the next section. 



6. Future High Resolution Observations 



Follow-up observations are important to check the predictions of our models. In the microlens- 
ing field, the idea of doing follow-up observations preceded the detection of planets, going back at 
le ast to 1998 in the ca se of MACHO-LMC-5, with the corresponding HST observations published 



m 



Alcock et al.l (120011 ). Then, it was successfully applied to derive more accurate parameters of 



the planetary systems, using HST or adaptive optics s ystem s at Keck or VLT, f o r instance in 



Bennett et al.l (|200J); bong et all (120091 ): Ijanczak et al.l (|2010h : Batista et al.l (J201l|); iKubas et al 



(|2012l ). In the future , it could be used to stud y planets in the h abitable zones of nearby dwarf 
stars, as suggested bv ldi Stefano fc Nightl (j2008h : Idi Stefanol 120121 ) . 



While the Bayesian analysis strongly favors a lens close to 1 M@, i.e, the upper limit permitted 
by the measurements, there is no airtight evidence against lower-mass lenses. Fortunately, this 
question can almost certainly be resolved by high-resolution imaging. The VIH source magnitudes 
are V a = 19.98 ± 0.03, I s = 18.71 ± 0.03, and H s = 17.35 ± 0.03. These come from CTIO H- 
band measurements, taken simultaneously with I and y-band observations. A linear regression 
between both fluxes gives a very accurate instrumental color of (I — H )ctio = 1-707 ± 0.002. After 
converting H to the 2MASS photometric system and / to the OGLE-III system by use of common 
stars in the CTIO field, the color becomes I — H = 1.36 ± 0.03. Subtracting this color from the I 
magnitude of the source in the OGLE-III system (18.71) returns the above value for H. 



-18- 



As we now show, the lens brightness must be similar to the source brightness. As can be 
judged from Figure O the field is relatively sparse, at a Galactic latitude of b = —5.095°, so these 
two stars (for the moment superposed) will very likely be the only two stars in their immediate 
high-resolution neighborhood. 

Let us consider three examples consistent with Equation ([9]), (M/Mq, -D^/kpc) = (0.1,0.4), 
(0.5,1.7), (1.0,2 . 8). T hese would have lens absolute magnitudes Mh = (9.6, 6.0, 3.4) according to 



Kroupa &: Tout! (|1997l ) and so H-^^ = (17.6, 17.1, 15.6). Of course, the extinction would be different 
at these different distances, but the entire column to the source is only Ah ~ 0.3. Therefore, a 
0.5 Mq lens star will be as bright as the source, and even a lens star at the bottom of the main 
sequence will produce an easily detectable amount of light (0.5 mag) over the expected source 
magnitude. 

Since all M/Dl combinations produce similar H magnitudes, such a measurement would, by 
itself, have little predictive power. But these various scenarios would yield substantially different 
J — H colors, which would add discriminatory power. Time has been allocated on various large 
telescopes to observe the field of this event, detect and measure the light coming from the source 
and lens stars. 

In addition, because of the relatively large proper motion, /i = 10.3 ± 0.8 mas/yr the lens and 
source could be separately resolved within about 5 years. This would then yield the angle of pro per 



motion (= tan 1 (vte i e/tte,n)) and so (from Figure [6|) the amplitude of tte (jGhosh et al. II2004I ). 



We can also consider follow-up observations with the Hubble Spa ce Telescope, which will be 



able to detect the lens-source relative proper motion as early as 2012 (jBennett et all 120071 ). The 
implied absolute magnitudes for the three examples given above (M = 0.1, 0.5, and 1.0 M©) are 
Mj = (12.3, 7.9,4.1), and the implied extinction free magnitudes are /l,o = (20.3, 19.1, 16.3). Since 
the extinction in the foreground of the lens is Aj < 1.0, this implies that the host star should be 
detectable with at least 9% of the /-band flux of the source star over the full range of main sequence 
host star masses. 

Finally, the ESA satellite GAIA to be launched in 2013 will image the Galactic Bulge down 
to V ~ 20 mag, so it may detect the lensing object and will certainly measure its proper motion if 
it does. 

In conclusion, these additional observations should reveal whether our choice of second-order 
effects (microlensing parallax and planet orbital motion) correspond to the reality. If we find 
contradictory results to our predictions, then it will be a strong argument to conduct the modeling 
of additional effects mentioned in Section 14.41 such as xallarap effect from a binary source, or 
involvement of a third body. 
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7. Final results, conclusions and perspectives 



In order to select the best among the three competing models (standard, parallax only, orbital 
motion and parallax), a simple comparison of x 2 values is not enough, because more refined models 
use more parameters. Taking into account these additional parameters by normalizing the x 2 
estimate by the number of degree of freedom is not the proper way to select the best model. A 
vast literature exists a b out m odel selection, and an application of different criteria to astrophysics 
is described in iLiddld (|20071 ). A simple way to take care of the larger number of parameters is 
to use the Akaike Information Criterion (AIC) (I Akaikd 1 1 974 ) . which introduces a penalty to the 



X 2 by adding twice the num ber of addition al parameters. Different criteria, such as the Bayesian 



Information Criterion (BI C) (|Schwarzlll978l ) or the Deviation Information Criterion, introduced by 
Spiegelhalter et al.l (|2002l ). can also be used. Formulas are given below, where Ak is the number 
of additional parameters, N the number of data points, and x 2 {@) the x 2 OI the average parameter 
set 9. 



AAIC = Ax 2 + 2 A A; 
ABIC = A x 2 + Ak log A 

ADIC = A(2^- x 2 (£)) (12) 



As there are 7 parameters in the standard model, 9 in the parallax only model and 11 in the 
orbital motion and parallax model, we see that our observed difference in x 2 01 8.6 in the parallax 
only model is only marginally significant according to the AIC (the expected difference is 4), while 
the observed difference of 54.4 (for the uq > 0) or 50.7 (for the uq < 0) is clearly an improvement 
of the orbital motion and parallax model over the standard one (the expected difference is 8). 
Similar results are obtained using DIC, while the difference between models returned by BIC is less 
significant. 

However, it is important to note that Ak strictly corresponds to the number of additional 
parameters only in the case of a linear regression problem. Here, we clearly have non-linear fits, so 
we should compute an "effective" number of parameters, which is difficult to estimate. The above 
conclusion should therefore not be taken as a quantitative one. 

A confirmation of the detection of second-order effects comes from the fact that in the orbital 
motion and parallax model, the degeneracy between uq < and Uq > models is clearly broken. It 
is instructive to plot both second-order effects vs each other, separately for the uq < and no > 
solutions. This is done in Figured where the orbital motion 7 is plotted vs. the parallax effect 7Te- 

If we remember that the ratio of the projected kinetic to potential energy must be smaller 
than 1 to get bound orbits, and that this ratio is proportional to 7 2 , where the proportionality 
constant depends on tte as given by Equation IE3|, it is easy to interpret these diagrams. When tte 
increases, the proportionality factor decreases, so that if 7 remains small enough, the bound orbit 
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condition is respected. In the uo > diagram (left side), small values of 7Te (say below 0.17) have 
a proportionality factor larger than 1. As 7 is nearly 1 for these solutions, they are ruled out, and 
the light curve confirms it, as few small x 2 solutions (black and red points) lie there. Larger tte 
are also excluded, as they correspond to large 7 values, although the light curve would favor such 
solutions. The only surviving region in this diagram is around 7Te ~ 0.35, corresponding to lens 
masses of half a solar mass, where the proportionality factor is about 0.4 and 7 slightly exceeds 1. 

In the uo < diagram (right side), although this is slightly disfavored by the light curve 
(Ax 2 = 6.7 for the best chain without the circular orbit constraint), there is a region where 7 is 
about constant at 0.7 for 7Te varying from to 0.4. The bound solutions correspond to the larger 
values of the 7Te domain (smaller proportionality factor), and they therefore agree with the range 
found in the uq > diagram. 

We therefore conclude that both solutions agree, and give bound orbits when the lens mass is 
about half solar, corresponding to a lens distance of about 1.6 kpc. 

If we now move to the post-bayesian analysis, we see that this solution favored by the light 
curve has some tension with the Galactic model constraint, because nearby lenses are rarer than 
more distant ones. But if we move to more distant lenses, we get many chains with unbound orbits 
or high eccentricities. By the way, the uq > solution, which has a lower \ 2 than the alternate 
uq < solution, is also the one where more chains correspond to unbound orbits. 

There is therefore a tension between Galactic and Keplerian priors, and the issue will only be 
solved photometrically, by measuring the light coming from the source and the lens. This will be 
the subject of a forthcoming article about this event. 

We conclude by giving the 1-D distributions of lens mass, lens distance, and planet orbit semi- 
major axis. The mass function for the lenses involved in these plots include main-sequence stars, 
brown dwarfs, but also white dwarfs, neutron stars and black holes, which may have large masses 
without violating the lens flux limit constraint. For the MS and BD stars, we adopt the following 
slopes of the present-day mass function dN / d\n{M): —0.3 between 0.03 and 0.7 M@, —1.0 between 
0.7 and 1.0 M , and —4.0 above. For the remnants, we adopt gaussian distributions, whose mean 
value, standard deviation, and fraction of total mass with respect to MS and BD stars below 1.0 M & 
are given in Table EJ 



For details about the choice of these numbers, please refer to iGouldl (|2000al ). 



In each diagram (see Figure [9] for mass and flux, and Figure [10] for distance and semi-major 
axis), the black curves show the full mass function, while the red curves show the mass function 
truncated at 1 M®. For MS stars, this limit is imposed by the lens flux constraint, and will be 
refined once we obtain the adaptive optics photometry of the individual stars in the field. WD at 
this mass are extremely rare; Jovian planets around pulsars (NS) have not been found, despite very 
extensive searches; and super-Jupiter planets orbiting BH are a priori unlikely. 



Let us first consider the lens mass distribution: the no-flux-limit (black) curve shows a huge 
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spike at expected NS position and a smaller bump corresponding to BH. Note that for these bumps, 
the uq < solution dominates, despite its x 2 handicap. This is because the Galactic model very 
strongly favors distant lenses, primarily because of the volume factor, and this overwhelms the 
modest preference of the light curve for nearby lenses. Because #e is roughly fixed, these distant 
lenses are massive. This preference is much stronger in the uq < solution, which can be seen in 
its rapid rise beginning at logM = —0.3. Note that the WD peak (at logM = —0.22) is clearly 
visible, especially in the uq > solution. 

The lens distance distribution basically looks at this same situation from the standpoint of 
distance. The new notable feature is that both MS and NS peaks are in the disk, while the BH 
bump is in the Bulge. And the semi-major axis distribution peaks at about 2-3 AU. 

^From these diagrams, we can estimate a most probable value of lens mass, distance and semi- 
major axis, and an asymmetric standard deviation read at 50% of the distribution corresponding to 
the red bold curves. We get a star and planet mass of M* = 0.671q'^3 M and m p = 1.5^ 3 Mjup, 
respectively, at a distance of D = 2.3 ± 0.6 kpc, and with a semi-major axis of a = 2±f AU. 

As a final note, it could be said that more complex models are worth exploring: the geometry of 
the caustic crossing, where the source passes close to the three-cusps tail of the caustic, is extremely 
sensitive to a third body (second planet or binar y companion to the l ens star). A similar g eometry 



where two planets were detected is described in iGaudi et al.l (120081 ) ; iBennett et al.1 (J201CJ) . These 



models could be investigated in a forthcoming paper, once we get the lens flux measurement from 
adaptive optics. 
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Table 2. Mass distribution of remnant stars (white dwarfs, neutron stars and black holes), with 
respect to main sequence stars and brown dwarfs below 1.0 M & 



Remnant < M > a ratio 



WD 

NS 

BH 



0.6 

1.35 

5.0 



0.07 22/69 
0.04 6/69 
1.0 3/69 
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Observations at Dome C 



ASTEP 400 is a 40cm Newtonian telescope installed at the Concordia base, located on the 

I). Although the aim of the project concerns 



Dome C plateau in Antarctica (D aban et al 



2011 



transiting planets, the ab ility to observe nea r-continuously during the antarctic winter and the 
excellent weather on site (jCrouzet et al.ll2010l ) imply that the telescope can usefully complement 
microlensing observations from other sites, even though the declination of the fields and their 
crowding make the analysis difficult. 

The observations with ASTEP 400 started on August 12, 2010, 15:27 UT after the first alert 
was sent by email and phone to Concordia. On August 14, 23:59 UT, the observations were stopped 
because the magnification had become too small for useful observations. The weather conditions 
were excellent. However, the seeing conditions were poor ( 3-4"), mostly due to the low-declination 
of the field and location of the telescope on the ground. The 229MB of data corresponding to 
500x500 cropped images were transmitted to Nice through satellite connexion around August 17- 
18 for an in-depth analysis. 

Although the data of this run were not good enough for being used in this study, this pioneering 
test will serve for improving the thermics of the acquisition system and get higher quality future 
observations of microlensing targets. 



B. Difference Image Analysis using pySIS 



The pySIS3.0 difference image anal ysis package is fully described in lAlbrow et al.1 ([2009). It 
is based on the orig inal ISIS package (|Alard fc Luptonl Il99«l : Ulardl hoooh but the kernel used 
to transform the reference image to the current image before subtraction is no longer analytic, 
but numerical. This allows dealing with images whose PSF cannot be assimilated to a sum of 
gaussian profiles. The numerical kernel has been introduced in image subtraction by IBramich 
(|2008l ). One of the regular problems encountered in difference image analysis is the choice of the 
best possible reference image. For astrometry, the best seeing image is generally a good choice, 
if the sky background is not too high. For photometry, it is our experience that stacking good 
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images improves the result, but only if these images have been acquired during a short time slot, 
to avoid light variations of the target and slow variations due for instance to small changes in the 



flat-fi eld. In order to choose good reference images, we use a suite of Astromatic software (jBertin 



namely SExtractor and PSFex. SExtractor builds catalogues of sources from all images of a 
given telescope, with their characteristics, and PSFex derives a model of the PSF of these images, 
from which we extract a few numbers to estimate the image quality (seeing, ellipticity, number of 
stars). This allows an almost automatic selection of the templates, with a final verification by eye 
to check the selected images. 

Once the images have been subtracted using these templates, the photometry of the target is 
done by iteratively centering the PSF on the light maximum close to the center of the image. This 
method enables the photometering of faint stars in the glare of nearby much brighter stars. 



C. Noise Model and Error Rescaling 

The standard procedure of rescaling error bars so that the x 2 /dof of each telescope data set is of 
order unity is acceptable if the resulting error bars roughly correspond to the dispersion of the data 
points at a given time for this data set. There is therefore an interplay between finding the correct 
model and rescaling the error bars, because a too large rescaling factor reduces the constraint from 
a given data set and allows the model to shift from the correct one. Our procedure has been to 
use rescaling factors which look plausible given the telescope size and site quality, trying to get 
a x 2 /dof of order unity only if the dispersion in successive data points is well reproduced by this 
rescaling scheme. This generally involves two parameters to modify the original photometric error 
bars e or i. One is a minimal error 6 m i n to reproduce the dispersion of very bright (or in this case, 
highly magnified) sources due to fundamental limitations of the photometry, such as flat fielding 
errors. The other is the classical multiplicative rescaling factor /. The adopted formula is: 

e = / V e ori + din (CI) 

The balance between both rescaling factors is given by comparing the cumulative distribution 
of x 2 /dof, ordered from magnification given by the model, to a standard cumulative distribution 
for gaussian errors. An example for OGLE and MOA distributions is given in Figure fTTT Table [3] 
gives the adopted values of both parameters for each telescope. It must also be noted that for 
some amateur telescopes, data were binned before this rescaling process. Finally, after the initial 
modeling was conducted, it was realized that three data sets could not accommodate the condition 
X 2 /dof ~ 1 (Auckland and VLO) or a positive source flux (Perth). They were therefore removed 
and the final models are based on 16 data sets. 
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D. Detailed Treatment of Limb-darkening Corrections 



It is generally difficult to treat the limb-darkening effect accurately, because it varies from one 
telescope data set to the other, primarily due to the different photometric bands involved. Second- 
order effects include hardware-specific variations of the spectral response (filters, CCDs), and for 
very broad bands, the atmospheric and interstellar extinctions. We adopt the linear limb-dark ening 
approximation using the formalism given by Equation EJ introduced by lAlbrow et al.1 (| 19991 ) and 
An et al.l (|2002i ). This is a common formalism in the microlensing community, but the more widely 
used formalism is based on the following equation 



/ A (cos0) =/ A (l) x [l-u A (l 



COS ( 



(Dl) 



where is the specific intensity at disk center and u\ is the linear limb-darkening coefficient. 

We compute the values of u\ from a stellar atmospher e mod el from the Kurucz ATLAS9 grid 



(Kurucz 



19931 ) using the method described by iHevrovskvl (|2007l ). As we converted the u\ values to 



our LDC parameter T\, we give the u\ values in Table EJ and recall the conversion relation 



2«A 
i - U\ 



(D2) 



For the source characteristics, we adopt the spectroscopic result (T e g = 6000 K, log g = 4.0, 
[Fe/H]=0.0). We thus get an estimate of the limb-darkening effect for each data set. Table [3] 
shows the adopted values of u\ for each individual telescope and band combination, both with 
and without interstellar extinction correction, together with the final number of data points after 
rejection of possible outliers, and details about the image reduction process. As can be seen, these 
coefficients vary even if the photometric band is supposedly the same. See for instance the Cousins 
/-band and SDSS i'-band, or the Wratten #12 filter (W12) used at some amateur telescopes to 
m imic an R band. More details about the derivation of these va lues for each telescope can be found 



m 



Hevrovskvl (|2007r ) : iFouque et all ll201(j ): Muraki et al.1 (|201lh . 



We have also computed the effect of interstellar extinction on the LDCs, as it is not a priori 
negligible: we typically find differences of a few 10 -3 for a data set with a filter and a few 10~ 2 
for unfiltered ones. The magnitude of the effect could then be neglected for filtered data sets, but 
not for unfiltered ones or very broad band filters. To give an idea, we compare its effect to the 
uncertainty of our spectroscopic determination of the effective temperature of the source, about 
150 K. For the OGLE /-band, the difference between LDCs with and without interstellar extinction 
correction for an adopted extinction of 1.2 mag corresponds to a shift in effective temperature of 
40 K, while for the Possum broad W12 filter, it corresponds to about 200 K, and for an unfiltered 
data set, it gives a shift of about 500 K. 
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Table 3. Data description for each telescope 



Observatory 


Band 


u\ (uncor.J 


u\ (ext. cor.) 


Binning 


AT 


f 




Photometry 




j 

1 










9 A 


n ni n 


DHTF DTA 

\ J v i 1 j 1 j 171 il 




j 

1 








-LOO 




n m n 

U.U1U 




CTIO 


V 


0.5931 


0.5908 




21 


2.0 


0.015 


DoPhot 


Auckland a 


W12 


0.5016 


0.4843 


Y 


5 


1 


0.003 


DoPhot 


FCO 


unfiltered 


0.5550 


0.5150 


Y 


4 


1.7 


0.003 


DoPhot 


FTS 


a 


0.4547 


0.4524 




68 


3.0 


0.002 


pySIS 


Kumeu 


W12 


0.5016 


0.4852 


Y 


17 


3.5 


0.003 


DoPhot 


Perth a 


I 


0.4325 


0.4297 




26 


3.0 


0.003 


pySIS 


FTN 


*' 


0.4547 


0.4524 




72 


3.8 


0.006 


DanDIA 


Possum 


W12 


0.5204 


0.5031 




77 


5.4 


0.003 


pySIS 


SAAO 


I 


0.4261 


0.4217 




2364 


2.2 


0.003 


pySIS 


SAAO 


V 


0.6114 


0.6086 




3 


1.2 


0.003 


pySIS 


VLO a 


unfiltered 


0.5442 


0.5022 


Y 


3 


1 


0.003 


DoPhot 


Wise 


unfiltered 


0.5522 


0.5131 


Y 


9 


0.8 


0.003 


DoPhot 


Canopus 


I 


0.4355 


0.4335 




50 


1.5 


0.010 


pySIS 


Danish 


I 


0.4394 


0.4370 




167 


2.2 


0.003 


pySIS 


MOA 


red 


0.4754 


0.4694 




3985 


1.05 





MOA DIA 


LT 


*' 


0.4543 


0.4513 




41 


4.7 





DanDIA 


Monet N 


I 


0.4414 


0.4392 




130 


4.3 


0.004 


pySIS 



Note. - 



a: This data set was not used in the final models 
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E. Detailed analysis of the Kepler constraint 



Figure [T2l shows several quantities plotted against 7Te,n- In each case the mean and standard 
deviation of all chain links within a given 7Te,n bin are calculated. The bottom panel shows the 
behavior of u) (= 7j_). Since 7Te,n is very similar to tte,±, this reflects the degeneracy between 
the perpendicular component s of Earth orbital m otio n (parallax) and lens orbital motion, which 
is analyzed in some detail by iBatista et al.1 (120111 ) and ISkowron et al.1 (120111 ). In the present case, 
the correlation is quite tight. The top panel shows the lens mass M = 6e/htte- Since 0e is nearly 
constant for different links in the chain, the mass scales M ~ tt^ 1 and so is peaked near 7Te 5 n = 
where 7Te is near its minimum (see Figured]). 

It is the two middle panels that enable one to understand why the orbital Jacobian favors 
relatively high masses. These show, respectively 711/70 and 7_i_/7o, where 



7|| 
7± 



s s, 
w, 

87T 2 / ^ s \ 3 
3 f E VT E + — I 



(El) 
(E2) 

(E3) 



Note that (7^ + 7y)/7o is the ratio of the so-called projected kinetic to projected potential 
energy, which has a strict upper limit of unity for bound orbits. We do enforce this limit, but the 
main impact of the orbital motion is more subtle. First look at 7_i_/7o- In the region away from 
7Te,n ~ 0, we have both 70 oc 7Te (Equation lE3|) and (very roughly) 7j_ oc tte (bottom panel). Hence, 
7_i_/7o is approximately constant in these two regimes. Much of the 7Ten > region is at or above 
the physical limit 7j_ =70, which is why these values are disfavored in Figure [6jc). 

Now examine T|| / To • By itself, 7y does not vary much with 7Te,n> so the form of the structure is 
basically just 7^ , which scales oc tt^ ^ away from zero. The point is, however, that the overall scale 
(which is set by the measurement of 711) is small, so that except near 7Te,n ~ 0, 7||/7o is extremely 
close to zero. Naively, this would seem to be disfavored by the virial theorem, but how does the 
Jacobian "know" about this? For simplicity of exposition let us consider circular orbits. For these, 
7||/7o = implies that the planet is exactly in the plane of the sky: either the orbit is exactly 
face-on (so this is always true), or the orbit just happens to be passing through the plane of the sky 
at the time of the event peak. The Jacobian is "unhappy" about either alternative because there 
is very little Kepler-parameter space relative to chain- variable space at such orbital configurations. 

Finally, note that in the immediate neighborhood of 7Te,n ~ (where M reaches its highest 
values), 7j_ is frequently at or above its physical limit (70), which is exacerbated by the relatively 
high values of 711 . This is reponsible for the modest suppression of extremely low parallaxes (high 
masses) in Figure El Thus, both Kepler and Galactic+fiux priors separately predict a relatively 
high lens mass, near the limit of what is permitted by the flux constraint. 
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Fig. 1. — Reduced spectrum with an overplotted synthetic spectrum for the adopted source char- 
acteristics, namely T e g = 5950 K, logg = 4.0 and solar metallicity. 
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Fig. 2. — Light curve of the microlensing event MOA 2010-BLG-477. Data from different observa- 
tories are distinguished by different colors. The upper panel shows the enlargement of the region of 
perturbation near the peak of the light curve. Note that the model light curve corresponds to the 
best-fit solution of the model considering both parallax and orbital motion and the corresponding 
parameters are presented in Table [TJ 
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Fig. 3. — Geometry of the lens system responsible for the microlensing event MOA 2010-BLG-477. 
In the lower panel, the filled dots represent the locations of the host star and planet. The closed 
figure composed of concave curves represents the caustic. The curve with an arrow represents 
the source trajectory with respect to the caustic. The big dashed circle centered at the host star 
represents the Einstein ring. The upper panel shows the enlargement of the region enclosed by 
a small box in the lower panel. The small filled circle represents the source star where its size is 
scaled by the caustic size. 
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Fig. 4. — OGLE-III 8x8 arcmin calibrated color-magnitude diagram (CMD) in I- and y-band. The 
red, yellow and green solid circles mark the centroid of the red giant clump, the blend, and the 
source position, respectively. 
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Fig. 5. — OGLE-III reference image of the BLG 176.8 field, in the region of our microlensed target. 
North is up and East to the left side. One pixel is 0.26". The red circle marks the position of the 
source (and the lens) from DIA analysis of CTIO images, in between two OGLE stars: the brighter 
of them (pointed at by a black arrow), at 0.37" from the source, corresponds in magnitude to the 
blended light measured by the model. The source is too faint to be detected in this crowded image. 
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Fig. 6. — Parallax vector (7te n, ^e.e) for the no > solution, from Markov Chain Monte Carlo fit 
to MOA-2010-BLG-477. Left panel (a) displays individual chain points, color-coded by A% 2 (1, 4, 
9, 16, 25, 36) for black, red, yellow, green, cyan, blue. Second panel (b) shows effect of weighting 
by the Galactic model, geometric Jacobian, and flux constraint. Lower parallaxes (and so higher 
masses) are heavily favored, except that too low parallaxes (corresponding to M > 1.0 Mq) are 
ruled out by the flux constraint. Third panel shows the effect of weighting using the Jacobian due 
to orbital motion parameters, which by itself disfavors the heaviest masses because these tend to 
imply unphysical orbits, but favors moderately high masses (see text). The final panel shows the 
effect of combining Galactic and Kepler priors. 
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Fig. 7. — Parallax vector (7Te,N) ^e.e) for the no < solution, from Markov Chain Monte Carlo fit 
to MOA-2010-BLG-477. See description of Figure [6] for details. 
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Fig. 8. — Variation of the orbital motion effect vs. the parallax effect for the orbital motion and 
parallax solutions. The left side corresponds to the no > solution, while the right side is for the 
uq < solution. 
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Fig. 9. — Distribution of lens mass and lens flux in ff-band (from left to right) computed from 
the post-bayesian analysis. In the left diagram, the black curves show the full mass function, while 
the red curves show the mass function truncated at 1 M & . The solid curves represent the no > 
solution, while the dashed curves correspond to the uq < one. The bold solid curves are the sum 
of the two other curves. In the right diagram, the curve corresponds to the red bold solid curve of 
the mass distribution. 
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Fig. 10. — From left to right: distribution of lens distance (in kpc) and semi-major axis of the 
planet orbit (in AU), computed from the post-bayesian analysis. For details, see the description of 
Figure [H 
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Fig. 11. — Cumulative x 2 /dof distribution for rescaled error bars in the case of OGLE photometry 
(left) and MOA photometry (right). 
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Fig. 12. — Binned lens characteristics as a function of MCMC variable 7Te.n- Bottom panel shows 
that 7Te,n is highly correlated with j±, the c omponent of orbital motion perpendicular to the 
projected planet-star found earlier by Batista et al.1 (120111 ). This is responsible for the 

poor determination of 7Te,n in Figure [6^a). Top panel: highest masses are associated with low 7Te,n- 
Third panel: these are disfavored because 7^/70 tends to saturate its physical limit (unity). Second 
panel: relatively high masses (with 7Te,n < 0) are favored because 711/70 (as well as 7_i_/7o) hover 
near 0.5, a value favored by the virial theorem. 



